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Abstract 

We propose the simplest possible renormalizable extension of the Standard 



O 

o 

r-| ■ Model — the addition of just one singlet scalar field — as a minimalist model 

for non-baryonic dark matter. Such a model is characterized by only three 

Oh! 

parameters in addition to those already appearing within the Standard Model: 



a dimensionless self-coupling and a mass for the new scalar, and a dimension- 
less coupling, A, to the Higgs field. If the singlet is the dark matter, these 
parameters are related to one another by the cosmological abundance con- 
straint, implying that the coupling of the singlet to the Higgs field is large, 
A ~ O(0.1 — 1). Since this parameter also controls couplings to ordinary mat- 
ter, we obtain predictions for the elastic cross section of the singlet with nuclei. 
The resulting scattering rates are close to current limits from both direct and 
indirect searches. The existence of the singlet also has implications for current 
Higgs searches, as it gives a large contribution to the invisible Higgs width for 
much of parameter space. These scalars can be strongly self-coupled in the 
cosmologically interesting sense recently proposed by Spergel and Steinhardt, 
but only for very low masses (< 1 GeV), which is possible only at the expense 
of some fine-tuning of parameters. 



1 Introduction 



It is an amazing fact of our times that even as our understanding of cosmology pro- 
gresses by leaps and bounds, we remain almost completely ignorant about the nature of 
most of the matter in the universe. According to recent fits to cosmological parameters 

dark matter of some sort makes up close to 30% of the total energy density. This 
is much more than what is inferred from inventories of the luminous matter we can 
see. Moreover, the successes of big bang nucleosynthesis suggest that only a fraction of 
this matter can be made of ordinary baryons, like massive compact objects, faint stars, 
etc.. Thus, unless gravity undergoes some drastic changes at distances larger than a 
few kpc (which is quite improbable from several points of view), we must postulate the 
existence of enormous amounts of dark matter of an unknown, non-baryonic origin. 

A simple argument shows why current theories of particle physics are so prolific in 
suggestions for what the nature of this dark matter might be ||. The vast majority of 
the proposed alternatives to the Standard Model involve new particles having masses 
which are of order m w ~ 100 GeV, and which couple with electroweak strength to the 
ordinary matter which we know and love. If any of these particles is stable enough to 
have a lifetime as long as the age of the universe, it makes a natural candidate for dark 
matter. It does so because its abundance is naturally of the right order of magnitude 
so long as its interaction cross sections have weak-interaction strength. The abundance 
comes out right because it is set by the annihilation rate for particles which are initially 
in thermal equilibrium with ordinary matter. Cosmologically interesting abundances 
follow pretty much automatically for particles whose mass is of order m w and whose 
annihilation cross sections have weak (or rather milliweak) interaction strength. (We 
give this argument in more detail within the body of the paper.) 

Better yet, supersymmetric models, which are perhaps the best motivated of the 
many theories which have been proposed, very often have such long-lived states, due 
to the natural existence there of a conserved quantum number, i?-parity, which keeps 
the lightest .R-odd state from decaying. These particles cry out for interpretation as 
dark-matter particles, and it is no surprise that these models are by far the most widely 
explored in the literature ||. 

Best of all, this explanation of the nature of dark matter can be tested experimen- 
tally. This is the direct goal of dedicated dark- matter detectors (4|, and an indirect 
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goal of accelerator searches for events with missing energy, showing that a weakly- 
interacting particle has escaped the detector. If Nature smiles on us we soon may be 
treated to the discovery of new physics in both of these kinds of experiments. In- 
deed, recently the DAMA collaboration has announced the detection of a dark matter 
signal, as indicated by their seeing an annual modulation of the counting rate in a 
Nal detector ||. However, the comparably precise data from the Ge detectors of the 
CDMS collaboration || do not support these findings. (These two experiments need 
not be in contradiction with each other if the spin-dependent part of the cross-section 
is enhanced relative to the spin- independent part @.) 

Our goal in the present paper is to present a slightly unorthodox view. Although 
very well motivated, supersymmetric models are very complicated and enjoy an enor- 
mous parameter space. This makes them unable to definitively predict what dark- 
matter detectors must see. Furthermore, unlike the extensive evidence for the exis- 
tence of dark matter, the arguments in favour of supersymmetry are almost exclusively 
theoretical. In our opinion, with the advent of good-quality data from dark-matter de- 
tectors, it behooves theorists to propose simple models for the dark matter which are 
consistent with present evidence, but which make definite predictions and so are easily 
falsifiable. These provide benchmarks against which other models and the data can 
be compared. We believe that it is only by comparing the implications of such models 
with one another, and with supersymmetry, that one can hope to properly interpret 
the data. 

The model we study in this paper was first introduced by Veltman and Yndurain 
P| in a different context. Its cosmology was later studied by Silveira and Zee M, and 



(with a complex scalar) by McDonald [[IIJ. It is the absolute minimal modification of 
the Standard Model which can explain the dark matter. It consists of the addition of 
a single spinless species of new particle, S, to those of the Standard Model, using only 
renormalizable interactions. To keep the new particle from interacting too strongly 
with ordinary matter, it is taken to be completely neutral under the Standard Model 
gauge group. Besides involving the fewest new states, the model is also just complicated 
enough to offer interestingly rich dark matter properties. Unlike the case if only spin- 
half or only spin-one singlet particles are added, it is possible for a singlet scalar to 
have both significant renormalizable self-interactions and renormalizable interactions 
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with some Standard Model fields. 

There is also a sense in which the model we propose is generic, should the dark 
matter consist of a single species of spinless particle. To this end, it is useful to ask the 
question of what a generic dark-matter model should look like. It is clear that the main 
property which one needs to ensure is the stability of the new particle, suggesting that 
the fields S i} representing these particles appears in the Lagrangian in even powers, so 
that its decay is forbidden. If this field Si is considerably lighter than the rest of the 
other exotic undiscovered particles, these may be integrated out, leaving an effective 
Lagrangian at electroweak scale which has the generic form 

C = kinetic term for $ + S*MSi + S*0 SM Si + ... (1.1) 

where the kinetic, mass terms and interactions with the SM (via the set of operators 
SM ) in general would depend on the spin of Si. The most important couplings at low 
energies are those of lowest dimension, corresponding to the lowest- dimension choices 
for the operators SM . Our model also has this form, with only a single singlet scalar 
S. In this language, our dropping of all nonrenormalizable interactions corresponds 
to keeping only those interactions which are consistent with (but do not require) all 
other exotic particles to be arbitrarily heavy compared with the weak scale. We might 
expect our model to therefore capture the physics of any more complicated theory 
whose impact on the dark matter problem is conveyed purely through the low-energy 
interactions of a single spinless particle. 

An additional, more tentative, incentive for formulating more models stems from 
recent indications of problems with subgalactic structure formation within the non- 



interacting cold-dark-matter scenario [pLT| . A 'generalized' form of cold dark matter 
may avoid these problems if its self-interactions^ can produce scattering cross sections 
of order 10~ 24 cm 2 in size ||12|| . Within the present context this proposal would require 
the masses of dark matter particles not to exceed 1 GeV. Within the minimal model 
described in this paper, we find this range of masses may be just barely possible, but 
requires unnatural fine tunings due to the relationship between masses and couplings 
imposed by the maintenance of the correct cosmic abundance of dark-matter scalars. 



1 The required self-interactions however also lead to spherical halo centers in clusters, which are 
inconsistent with the ellipsoidal centers indicated by strong gravitational lensing data fu| . 
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This paper is organized as follows. In the next section we identify the three pa- 
rameters which describe the model, and determine the general conditions which lead 
to acceptable masses and to sufficiently stable dark matter. In section 3 we calculate 
the annihilation cross section of S'-particles and give the resulting cosmic abundance as 
a function of masses and couplings. This calculation is similar to the analysis of Ref. 



10]. We perform the numerical analysis for the most interesting part of the parameter 



space, with 100 GeV < m h < 200 GeV and 10 GeV < m s < 100 GeV. In section 4 
we obtain the cross section for elastic scattering with ordinary matter and apply the 
constraints, imposed by direct and indirect searches. Section 5 computes the cross 
sections for the missing energy events which are predicted for colliders due to the pair 
production of S particles. It also contains a prediction for the degradation of the Higgs 
boson signal at hadronic colliders, when the Higgs boson is allowed to decay into a pair 
of 5* particles. Our conclusions are reserved for section 6. 



2 The Model's Lagrangian 

The lagrangian which describes our model has the following simple form: 

C = Csm + \d»Sd»S -!^S 2 -^S'-\S 2 H^H, (2.1) 

where H and C SM respectively denote the Standard Model Higgs doublet and la- 
grangian, and S is a real scalar field which does not transform under the Standard 
Model gauge group. (Lagrangians similar to this have been considered as models for 
strongly-interacting dark matter WM and as potential complications for Higgs searches 



15 1 . The same number of free parameters appears in the simplest Q-ball models fig] .) 
We assume S to be the only new degree of freedom relevant at the electroweak scale, 
permitting the neglect of nonrenormalizable couplings in eq. ( [2.i| ), which contains all 
possible renormalizable interactions consistent with the field content and the symmetry 

s -> -s. 

Within this framework the properties of the field 5* are described by three param- 
eters. Two of these, A s and are internal to the 5* sector, characterizing the S mass 
and the strength of its self-interactions. Of these, \ s is largely unconstrained and can 
be chosen arbitrarily. We need only assume it to be small enough to permit the pertur- 
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bative analysis which we present. Couplings to all Standard Model fields are controlled 
by the single parameter A. 

We now identify what constraints are implied for these couplings by general con- 
siderations like vacuum stability or from the requirement that the vacuum produce an 
acceptable symmetry-breaking pattern. These are most simply identified in unitary 
gauge, \f2H^ = (h,0) with real h, where the scalar potential takes the form: 

V = f S 2 + ^ S 2 h 2 + ^S* + ?f(h 2 - v 2 EW ) \ (2.2) 

\h and v EW = 246 GeV are the usual parameters of the Standard Model Higgs potential. 

1. The Existence of a Vacuum: This potential is bounded from below provided 
that the quartic couplings satisfy the following three conditions: 

A s , A A > and (2.3) 
A s \h > A 2 for negative A. 

We shall assume that these relations are satisified and study the minima of the scalar 
potential. 

2. Desirable Symmetry Breaking Pattern: We demand the minimum of V to have 
the following two properties: It must spontaneously break the electroweak gauge group, 
(h) 7^ 0; and it must not break the symmetry S — > —S, so (S) = 0. The first of these 
is an obvious requirement in order to have acceptable particle masses, while the second 
is necessary in order to ensure the longevity of S in a natural way. (S particles must 
survive the age of the universe in order to play their proposed present role as dark 
matter.) 

The configuration h ^ and S = is a stationary point of V if and only if v 2 EW > 0, 
in which case the extremum occurs at hl xt = v\ w . This is a local minimum if and only 
if 

m 2 Q + Xv 2 EW >0. (2.4) 

A second local minimum, with h cxt = and S 2 xt = —ml/X s , can also co-exist with the 
desired minimum if A > and A 2 < A/jA s . This second minimum is present so long 
as itLq < and — Am 2 , > X s XhV EW . Even in this case, the minimum at S cxt = and 
/i 2 xt = v EW is deeper, and so is the potential's global minimum, provided that 

< -m 2 < v 2 EwV ^X~ s . (2.5) 
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Throughout the rest of this paper, the above conditions are assumed to hold, so that 
the model is in a phase having potentially acceptable phenomenology. It is therefore 
convenient to shift h by its vacuum value, h — > h+v EW , so that h represents the physical 
Higgs having mass m 2 h = \hV 2 EW . The ^-dependent part of the scalar potential then 
takes its final form 



and the S mass is seen to be m 2 = m 2 , + \v 2 w . Our prejudice in what follows is that 
this mass lies in the range from a few to a few hundred GeV, in which case the resulting 
dark matter will be cold. 

3 Constraints from Cosmological Abundance 

We next sharpen the cosmological constraints on the model by demanding the present 
abundance of S particles to be close to today's preferred value of Q s h 2 . This imposes 
a strong relationship between the parameters A and m s , which we now derive. 

We start by assuming that the S particles are in thermal equilibrium with ordi- 
nary matter for temperatures of order m s and above. This is ensured so long as the 
coupling A is not too small. Just how small A must be is determined by the following 
argument. Thermal equilibrium requires the thermalization rate, T t h, to be larger than 
the universal Hubble expansion rate, H. The constraint on A comes from the demand 
that this be true throughout the thermal history of the universe, down to temperatures 
T < m s ~ m w . But T t h and H vary differently with time as the universe expands, 
because they differ in their temperature dependence. On one hand, in the radiation- 
dominated epoch which is of primary interest to us, H ~ T 2 /M P1 , where M P1 denotes 
the Planck mass. On the other hand, the thermalization rate varies as r th ~ A 2 T for 
T 3> m h , and as r th ~ A 2 T 5 m^ 4 for T <C m h . These temperature dependences imply 
the ratio R t h = ^th/H is maximized when T ~ rrih ~ m w , taking the maximum value 
-Rth ~ A 2 Mpi/m ft . S particles are therefore guaranteed to remain thermalized (or get 
thermalized) down to the electroweak epoch, T ~ m h ~ m w , if this maximum ratio is 
required to be of order one or larger, implying 




(2.6) 




(3.1) 
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Figure 1: The Feynman graph relevant to 5-particle annihilation via Higgs exchange. 
Various annihilation channels are open or forbidden, depending on the value of 2m s . 

Once thermalization is reached (and in the absence of S decays), as we shall hence- 
forth assume, the primordial S abundance is determined by the S particle mass and its 
annihilation cross section. This cross section depends very strongly on the unknown 
Higgs mass, and on which annihilation channels are kinematically open. 

An independent, interesting issue is the fate of a scalar condensate that might sur- 
vive from the inflationary epoch. After the Hubble rate drops below ms, coherent time 
oscillations of the singlet field begin. These oscillations can be regarded as the oscilla- 
tions of a Bose condensate of S particles which is not in thermal equilibrium with other 
matter. The fate of the condensate depends on the initial value of the S field and two 
possibilities must be distinguished. If the initial value of the condensate is sufficiently 
small so that the energy density in the oscillations, ~ m^(S) 2 , is smaller than the 
energy density of radiation, ~ T 4 , then the thermalization of this condensate occurs 
exponentially fast. The rate is given by A 2 T or A|T, whichever is larger. When the 
initial value of the condensate is of the order of the electroweak v.e.v., the condensate 
will therefore completely disappear if A or X s is larger then y/Mw/Mpi, just as in the 
thermalization condition (3.1). The situation is quite different when the initial value 
of the S field is very large (Mpi, for example). The S condensate then dominates the 
energy density in the Universe and it behaves exactly as the inflaton condensate. The 
absence of the direct decay of S particles in this case may prevent the universe from 
reheating [[17]] 0. In this paper we assume that S field does not drive inflation, and we 
limit ourselves to the first possibility. 

Since the temperature domain for which annihilation is most important is T ann ~ 
2 We thank Lev Kofman for pointing out this possibility 
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0.05 m s , it is the nonrelativistic annihilation cross section which is relevant. In our 
model the expression for the annihilation rate depends on which phase within which it 
occurs. If it occurs within the Higgs phase, i.e. if T ann is low enough so that it occurs 
after the electroweak phase transition, the result is given by evaluating the tree-level 
graph of fig. (1) for s-channel annihilation, SS — > X, which in the nonrelativistic limit 
gives 

_ „ _ 8A 2 u| w jp /oo\ 

a ^ Vrel ~ {Aml-mlY + mlTl^ ^ 



with F x := lim 



Here I\ is the total Higgs decay rate, and T- hX denotes the partial rate for the decay, 
h — > X, for a virtual Higgs, h, whose mass is m~ h = 2m s . Eq. (|3.2| ) also assumes that 
m s < rrih, so that direct (contact) annihilation to a pair of physical Higgses via the 
XS 2 h 2 interaction term is forbidden. 

Of particular interest are the large- and small-m<j limits. The small-m s limit of 



eq. (|3.2|) - m s m w , m h - implies the asymptotic behaviour: 

A 2 ml 



Oann^reZ OC if m s < m h , (3.3) 

m h 

and the coefficient of proportionality depends strongly on the accessibility of certain 
decay channels (X = tttt or and so on). 

For large m s the dominant contributions to the annihilation cross section come from 
W^VF - , ZZ and hh final states. (The latter originates from the XS 2 h 2 interaction 
term, whose contribution must be summed with ( |3.2| )). Neglecting terms which are 
0(v 2 EW /m 2 s ) in the result, we find the large-m s behavior of the annihilation cross section 
to be 

A 2 

OaxmVrel ~ |^2~ > if m s > m h- (3-4) 

These asymptotic forms are useful in what follows for understanding what the cos- 
mological abundance constraint implies for the coupling A in the limit where m s is 
very large or very small. Our results for the annihilation cross section agree with the 
calculation of ref. p0[] . 
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Figure 2: Four samples of the logA-m s relationship between A and m s , which gives 
the correct cosmic abundance of S scalars. For these plots the Higgs mass is chosen to 
be 100, 120, 140, and 200 GeV. The abundance is chosen to be Q s h 2 = 0.3 



These expressions may be used with standard results for the Standard Model Higgs 
decay widths to predict how the primordial S'-particle abundance depends on the pa- 
rameters m s and A. Standard procedures fl8| give the present density of S particles 
to be 

= d-°7-i°V, , (3 . 5) 

gl /2 M P1 GeV (av ie} ) 

Here g*, as usual, counts the degrees of freedom in equilibrium at annihilation, Xf — 
m s /Tf is the inverse freeze-out temperature in units of m s , which can be determined 
by solving the equation Xf ~ In 0.038(g*x/) -1 / 2 Mpi m s (av ) . Finally, (• • •) denotes 
the relevant thermal average. 

The abundance constraint is obtained by requiring Q s h 2 to be in the cosmologically 
preferred range, which imposes a relation between A and m s - For our numerical results 
we use Q s h 2 = 0.3, which corresponds to a large value for Q s = 0.6, perhaps the largest 
possible value consistent with observations. We choose such a large Q s in order to 
be conservative in our later predictions for the signals to be expected in dark matter 
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detectors, since larger f2 s corresponds to smaller A. Choosing instead the central value, 
Q s h 2 ~ 0.15, obtained from various cosmological fits would give values for A 2 which 
are about twice as large as those which we use in what follows. Since it happens that 
Xf ~ 20 is approximately constant for the range of cross-sections and masses expected 
in this model, the abundance condition is equivalent to holding (av) constant as A and 
m s are varied (provided is held constant). 

Fig. (2) plots the relationship between A and m s which is predicted in this way by 
the requirement that Q s h 2 ~ 0.3. For most values of m s this curve is well described by 
the above simple formulae, which give sufficient accuracy in most parts of the parameter 
space defined by varying m s , A and m^. Important exceptions to this statement apply 
in kinematically special regions, such as the Higgs threshold (2m s ~ rrih) and two- 
particle thresholds in the final states (2m s ~ 2m?, or 2m s ~ 2m w , and so on), where 
a more sophisticated treatment is required. In our quantitative work we follow the 



numerical procedure outlined in Ref. |19| in these special regions. 

For m s sufficiently large or small, the asymptotic expressions ( |3.3| ) and ( |3.4| ) show 
that the abundance constraint forces A to become large, eventually becoming too large 
to believe perturbative expressions like eq. G3.2| ). In particular, if annihilation should 
occur before the electroweak transition, then the asymptotic relation between A and 
m s becomes: 

A ~ (3.6) 
lOTeV V ' 

so demanding the perturbative regime (A < 1) gives the upper bound m s < 10 TeV. 

For small m s we consider the case m s = 500 MeV, for which the dominant anni- 
hilation channels are 7r + 7r~, 7r°7r° and P0|. In this case Q s h 2 < 0.3 is achieved 



if the coupling satisfies the constraint A > 2. Smaller values of A for this range of m s 
would lead to over-abundant scalar particles and an overclosure of the universe. 

There are several important points concerning the abundance condition which bear 
emphasis: 

1. For all m s < few TeV (and away from poles and particle thresholds) the abun- 
dance constraint requires A ~ O(0.1 — 1). In this sense this model of dark matter 
is "natural", in that obtaining the right primordial abundance does not require 
any fine tuning or special choice of the parameters. (This property is shared by 
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many models for which the dark matter is a weakly-interacting species of particle 
whose mass is of order m w .) 

2. The coupling A has to be significantly suppressed (down to the level of 10~ 4 — 10~ 3 ) 
near the Higgs pole. This is because the Higgs resonance is rather narrow, and 
this narrowness considerably enhances the S annihilation rate, especially if 2m s 
is slightly smaller than [19| . 

3. Different decay channels dominate the total annihilation cross section for different 
ranges of m s . However, the range of values of most experimental interest lies 
between the b and W thresholds, for which it is the bb final state that is most 
important. 

4. Since the abundance constraint is concerned with the strength of the interactions 
between S scalars and ordinary matter, it is largely independent of the strength 
of the S self-coupling, A s . This leaves A s completely free to be adjusted. Unfor- 
tunately, although the S particles therefore can be very strongly interacting, this 
in itself does not make them useful to solve the recently-perceived problems with 
galaxy formation [11|. This is because the solution of these problems requires 



interaction cross sections which are of order 10 24 cm 2 , and cross sections this 



large require m s < 1 GeV in addition to large A s [O, P|. Unfortunately masses 



this small require fine tuning in this model, due to the relation m| = m 2 , + \v 2 EW . 
Since, as we saw earlier, the abundance constraint requires A to be of order one 
or larger for small m 3 , we require a part-per- million cancellation between Xv EW 
and m 2 , in order to obtain small values for m 2 . 

It would be interesting to explore whether more natural possibilities are possible 
with a less draconianly minimal scalar model. One possibility might be enhance 
the cosmic abundance using couplings to SM fermions, while at the same time 
suppressing the coupling to the Higgs to avoid the problem of fine-tuning. An 
alternative approach is to arrange for relatively light, but strongly self-interacting 
scalars whose coupling to visible matter is suppressed below the level required by 
the thermalization argument which led to eq. ( |3.1| ), and instead arrange for the 
correct cosmic abundance of S'-particles in another way, perhaps as the result of 
an earlier inflationary epoch 
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4 Implications for Dark Matter Detectors 



The connection between A and m s derived from the abundance constraint in the pre- 
vious section is very predictive. In this section we explore its consequences for current 
dark matter searches, while the next section describes implications for Higgs physics 
in collider experiments. 

The sensitivity of dark matter detectors to S particles is controlled by their elastic 
scattering cross section with visible matter, and with nuclei in particular. This cross 
section enters in one of two ways: (z) it is the directly relevant quantity for experiments 
designed to measure the recoil signal of dark-matter collisions within detectors 0, |6], |22, 



23| . |24| . |25| ; (ii) it controls the abundance of dark matter particles which become trapped 
at the terrestrial or solar core, and whose presence is detected indirectly through the 
flux of energetic neutrinos which is produced by subsequent S-particle annihilation 



26, 27, 28, 2S 



N ' N 

6 

^hNN 



Figure 3: The Feynman graph relevant to scalar-nucleon elastic scattering. 

The Feynman diagram describing elastic S-particle collisions with nucleons and 
nuclei is given by t-channel Higgs exchange, as shown in fig. (3). If we write, for 
slowly-moving spin- J nuclei, the relevant nuclear matrix element as 



2 J + _ . 

spins / 



then the nonrelativistic elastic scattering cross section obtained by evaluating fig. (3) 
becomes 

„ J = ftWf^, (4 . 2) 
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where m* = m s m N / (m s + m N ) is the reduced mass for the collision. 

Modelling the nucleus as A independent nonrelativistic nucleons leads to the ex- 
pectation that ^.nucleus ~ ^Amcieon , making it convenient to relate the elastic cross 
sections for nuclei and nucleons by 

cr c i (nucleus) = - a c i (nucleon), (4.3) 

mi(p,S) 

The Higgs charge of the nucleon can be related to the nucleon mass and the trace 



anomaly, following ref. [ITJ|. The result is sensitive to the mass of the strange quark 



and its content in the nucleon in the + channel. Taking the strange quark mass to be 



170 MeV and (iV|ss|./V) ~ 0.7, (see, for example, Ref. [plfl ), we deduce the estimate 



_ 340 MeV 

^nucleon — QhNN ~ j l^- V 

Vew 



leading to 

del (nucleon) 



A 340 MeV A 2 / m 



p 



m\ ) \n(m p + m s ) 

100 GeVA 4 /50 GeVA 2 / in _ 42 2 \ , , 
j I J (^20 x 10 42 cm 2 j. (4.5) 



4.1 Recoil Experiments 

It is instructive to see what kind of cross section is implied by the primordial cosmic 
abundance relation, A(m s ). Using = 120 GeV and m s = 40 GeV, the abundance 
condition 0.1 < Q s h 2 < 0.3 implies 0.05 < A < 0.085. Using these in eqs. ( |4.5| ) then 
gives 5 x 10~ 44 cm 2 < <j e i(nucleon) ~ 1.5 x 10 -43 cm 2 , which is about a factor of 20 
smaller than the range of nucleon scattering cross sections which are accessible in the 
DAMA and CDMS experiments [|, ||. 



Fig. (^|) shows how the scattering cross section we obtain in this fashion is related 
to current experimental limits, under the standard assumption that the mass density of 
S'-particles in our galactic halo is 0.3 GeV/cm 3 , having velocities of order 200 km/sec. 
This plot shows the lowest possible values of the cross sections, as it uses the maximum 
possible abundance Q s h 2 = 0.3. For smaller values of the cosmic abundance, A 2 is 
higher and the cross sections could be up to three times larger than those shown in 
fig. (EI). The predictiveness of the A — m s relation makes this model much easier to 
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Figure 4: The predictions for the elastic cross section, function of m s , which 

follows from the X(m s ) dependence dictated by the cosmic abundance. Also shown by 
a dashed line is the exclusion limit from the CDMS experiment |J . 

falsify than are more complicated models, with much of the parameter space covered 
by the next generation of experiments Most importantly, the projected sensitivities 
of the CDMS-Soudan and Genius experiments will completely cover the range m s < 50 
GeV, for values of the Higgs mass between 110 and 140 GeV. As we show in the next 
section, this range of masses and coupling constants has important implications for the 
Higgs searches at colliders. On the other hand, there exists the possibility of completely 
"hiding" the dark matter by choosing OArrih < m s < 0.5m/j. In this case annihilation 
at freeze-out is very efficient, requiring small A's which lead to elastic cross sections 
suppressed to the level of 10~ 48 cm 2 . These levels of sensitivity to a c i(nucleon) are not 
likely to be achieved in the foreseeable future. 

Our model of a singlet real scalar predicts a smaller signal for underground detectors 
than does a model where the dark matter consists of iV singlet scalars (including the 
model considered in ref. jl(|, for which N = 2). This is because the abundance of every 
individual species must be 1/N of the total dark matter abundance, fli = Q tot /N. This 
requires a larger annihilation rate at freeze-out for every species, and so an enhancement 
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of the coupling A by y/N, relative to the values of coupling constant calculated in this 
paper. By contrast, since the local halo density does not depend on N, the signal from 
an iV-component scalar dark matter model will be N times larger than what is found 
here for the signal from our single-component dark matter model. 



4.2 High-Energy Neutrino Searches 

The elastic cross section calculated above, a e \, also controls the expected flux of high- 
energy neutrinos which would be emitted by S particles which are captured at the 
centre of the sun or the earth. It does so because of the following scenario, which 
describes how the abundance of such particles is determined. 

5* particles which lose enough energy through scattering get trapped in the gravita- 
tional field of the sun or the earth. Further scattering permits them to further dissipate 
their energy, until they eventually accumulate at the solar (or terrestrial) center. Their 
density at the centre grows until the accumulation rate precisely balances their rate of 
removal due to pair annihilation, leading to an equilibrium abundance. Because of this 
balance the total rate of annihilations may be computed given the capture rate, and 
so also given the elastic scattering rate. 

The detection of these captured particles is based on observing the high-energy 
neutrinos which are among the production products of these annihilations. These 
neutrinos can escape from the solar (or terrestrial) centre, and can be detected by 
neutrino telescopes, which look for the energetic muons which are produced when the 
neutrinos interact with rock or ice in the detectors' immediate vicinity. In this way a 
signal is predicted for detectors like Baksan, Kamiokande, Macro, Amanda and others. 



To predict this signal we use the results of Ref. [33], which follows the original 
treatment in Refs. [33], and estimate the neutrino flux at the surface of the earth due 
to S-particle annihilation in the sun as 

4> vQ ~ (560CHTV 1 ) iV cff (7 el) 36 (4.6) 

m s 

In this formula N e g is the average number of neutrinos produced per annihilation event, 
and o"oi,36 is the S-proton elastic scattering cross section in units of 10~ 36 cm 2 . 

The determination of N e s requires a study of the final states which are available in 
S annihilations, which are mediated by s-channel Higgs-exchange, as in Fig. (1). The 
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number of neutrinos produced depends on the value of m s , which controls whether the 
main annihilation products are W + W~, bb or other light hadrons. (Direct branching 
to a pair of neutrinos is obviously very small for Higgs- mediated decays.) The total 
production of energetic neutrinos turns out to be quite significant, as all decay products 
typically produce final state neutrinos (and so non-vanishing N e g) due to cascades of 
weak decays. The resulting value for iV eff and the energy spectra of the produced 
neutrinos have been meticulously simulated by Ritz and Seckel in Ref. [ JJ| and more 
recently by Edsjo pEfl . 

To obtain an upward-going flux of energetic muons, one must compute the proba- 
bility that a neutrino directed towards the detector produces a muon at the detector. 
This probability, P(E), is a quadratically rising function of neutrino energy, E, which 
must be convoluted with the computed neutrino fluxes [57, [3j|. Approximating this 
probability as P(E) ~ l(r 8 (£/100 GeV) 2 , we use the results of Ref. for 60 GeV b 
quarks injected into the center of the sun, which predicts N e Q to be 0.3. The second mo- 
ment of the neutrino energy distribution, N e g ((E/m s ) 2 ) becomes 0.006. Extrapolating 
these results to the case m s = 50 GeV gives 

i _ 14 ( cr c i(nucleon) ^ __ 2 „-i 
10~ 44 cm 2 



°- 8 x 10 -rrbz — t- cm s ■ ( 4 - 7 ) 



This should be compared with experimental limits on the flux of energetic upward 
muons (as obtained, for example, by the Kamiokande, MACRO and Baksan experi- 
ments, [|2B], |27J, ^3]), which is 4>n ^ 1-4 x 10~ 14 cm~ 2 s _1 . Using the value of <7 e i(nucleon) 
suggested by abundance calculations (c.f. eq. ( |4.5|) ), we can see that the predicted flux 
is right at the level of current experimental sensitivity. The constraints on the flux 



coming from the center of the earth can be equally important [10 



There are, of course, caveats to the blanket use of these constraints, since loopholes 
can exist for some values of the parameters. For instance, according to Ref. [ 35[| , the 
average energy of the neutrinos is about 10% of m s , and so for lower values of m s 
the muons passing through the detectors might be close to or below the experimental 
thresholds. We conclude that although the prospects are good for observable signals 
in these detectors, more detector-oriented studies need to be done in order to exploit 
fully the limits which may result. 
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5 Implications for Collider Experiments 



We next turn to the implications of the model for Higgs searches at colliders. As we 
saw from the primordial abundance constraints, over most of parameter space A ~ 
O(0. 1) — 0(1). This means that real or virtual Higgs production might frequently also 
be associated with S production, potentially leading to strong missing energy signals. 



(See refs. [15] for discussions of the implications of related scalar models for collider 



experiments and refs. jHJ for the discussion of invisible Higgs decay in supersymmetric 
models.) Since the extent to which this actually happens depends strongly on the 
S'-particle mass, we consider the main alternatives successively. 

5.1 2m s < nih 

This is the most interesting possibility, since it permits the decay of the Higgs into a 
pair of S'-particles. The decay width, calculated at tree level, is given by 



A 2 ?; 2 / 4m 2 

T h ^ ss = ^^Jl-^f for 2m s <m h . (5.1) 
8nm h y m z h 

Needless to say, since any produced S particles are unlikely to interact within the 
detectors, the large values for A suggested by cosmic abundance can give a very large 
invisible decay width to the Higgs, without changing the Higgs' other Standard Model 
couplings. This fact obviously has many consequences for Higgs searches at colliders, 
since Higgses will be produced at the rates expected in the Standard Model, but will 
mostly decay invisibly into S pairs. To quantify this observation we calculate the ratio 
of the above invisible width to the width for the decay h — > t + t~: 



T h -> ss \ 2 v% w I Ami 



T h ^ T+T - m\m 2 T Y m\ 

Using the usual cosmic abundance constraints, obtained in Section 3, one can see that 
this ratio can be as large as 1000, so that the decay of the Higgs into a pair of S 
scalars is by far the more probable event. More useful information is given, however, 
by another ratio, 

Br h ^ T+T -{SM + S) T 

h, total 

(SM) 

Br^ T+T - (SM) ~ T h ^ ss + f\ total (SM) ' l " J 
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Figure 5: The ratio, R, of the total Higgs width in the Standard model over the same 
width in the Standard Model supplemented by the singlet scalar, plotted as a function 
of m s . 

where the last equality uses the fact that Higgs production rates are not affected by 
the existence of a new scalar. R quantifies the deterioration (relative to the SM result) 
of the expected signal for Higgs decaying into visible modes due to the adding of the 
new scalar. We plot R, in fig. (5), against m s up to the Higgs mass for the same values 
of rrih and A as in fig. (2). 

As is clear from the plot, for m^ = 100, 120 and 140 GeV and 2m s < m^ the 
invisible width dominates the total width everywhere except in the immediate vicinity 
of 2m s = rrih. R shrinks near this point for two reasons. First, m^ near 2m s is close to 
threshold for producing two S particles in h decay, and so the invisible rate is phase- 
space suppressed in this region. More importantly, the size of the coupling, A, allowed 
by abundance arguments is smaller for m s in this region, due to the enhancement of 
the primordial annihilation cross section due to proximity with the Higgs pole. 

However, the plot also shows that the presence of the invisible width already down- 
grades the expected signal by a factor of 10 or more when m$ — 0.3m/,, GeV for 
mh = 120 or 140 GeV. This means that over 10 times more Higgs particles must 
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be produced in order to reach the same level of signal in all visible decay modes, 
as compared to what is expected purely within the Standard Model. This implies a 
tremendous suppression of the observable Higgs signal at the LHC and especially at 
the Tevatron — possibly even precluding its discovery at these machines (see, however, 
ref. [41], El' HI f° r a more optimistic point of view). The main Higgs production re- 



action in this case would be two jets plus missing energy, a process having enormous 
backgrounds in hadron machines. 

On the other hand, for rrih close to or above the W + W~ threshold, the similarity 
in size between A and the electromagnetic coupling, e, implies that the invisible width 
does not completely dominate other decay modes. In such a case (fig. (5), m h =200 
GeV) the decline in R is not dramatic, and the existence of an invisible signal need not 
preclude finding the Higgs in a hadron machine. 

More can be said about an invisibly decaying Higgs at e + e~ machines. Although the 
missing energy signal would be missed in an orthodox Standard-Model Higgs search, 
such as one based on 6-tagging, more model-independent searches have been performed 
for scalars produced with Standard Model cross sections, but decaying invisibly. These 
searches at LEP exclude such an invisibly-decaying Higgs unless rrih > 106.7 GeV [13 . 



5.2 2m s > nrih 

For this mass pattern, S particles cannot be produced by real Higgs decays, and so 
arise only through virtual Higgs exchange. Again, once produced, these particles are 
not expected to interact within the detector and so look like missing energy above an 
energy threshold, E > 2m s . 

A missing-energy signal can be searched for at e + e~ colliders, where the back- 
grounds are very well understood. Since hadronic machines are unlikely places for 
seeing such a signal, we do not consider them further here. 

Fig. (6) shows the main mechanism for the production of two S scalars at LEP (or 
at a linear collider). The differential cross-section for this process which results from 
the evaluation of this graph is as follows: 

s + (t — m 2 z )(u — m 2 z )/m 2 z 



da 



dt du 



16(2 



V I 



SnnOq, 
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Figure 6: The dominant mechanism for pair-producing S'-particles at LEP, via s- 
channel Higgs exchange. 



x 



s + t + u 



mi 



mi 



4m| 



(s + t + u 



mi. 



(5.4) 



Here we follow the notation of ref . |39| , with g L = — \ + s 2 w and g R = parameterizing 
charged-lepton couplings to the Z boson, a z = em z / (c w s w ) being the effective coupling 
of the Higgs to Z boson, and s w and c w denoting the sine and cosine of the Weinberg 
angle. The Mandelstam variables, s, t and u, are defined as for a 2-body to 2-body 



process, with their sum given by s + t + u 



m: 



+ Q 2 , where m z is the Z-boson mass 



and Q 2 is the square of the invariant mass of the two S'-particles. 

Expression ( |5.4|) also applies to the previous case, where 2m s < m^, in which 
case the cross section is dominated by the Higgs pole, corresponding to real Higgs 
production. It is clear that the constraint on the missing energy implied by eq. (|5.4| ), 
and the experimentally measured values for such a process at LEP, cannot lead to a 
limit on the invariant mass of two S particles better than 100 GeV for a Higgs mass 
of around 100 GeV. For larger values of m^ a possible bound on the invariant mass of 
the S pair is considerably relaxed. 



6 Conclusions 

We have presented the first study of the minimal model for non-baryonic dark matter, 
which consists of the Standard model plus a singlet scalar. This model is characterized 
by one additional real scalar field and three new parameters. The absence of linear 
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and cubic terms in S is required to ensure that the the new scalar is sufficiently stable 
to contribute significantly to the dark matter currently present in the Universe. This 
stability also precludes the development of a nonzero v.e.v. for S. The study of the 
model's scalar potential shows that electroweak symmetry is spontaneously broken 
while (S) vanishes for a significant domain of parameter space. 

At the renormalizable level the three parameters describe the S particle mass, its 
self-coupling and its couplings to other Standard Model fields, which all are mediated 
by a coupling to the Standard Model Higgs, of the form XS 2 (H^H). The simplicity 
of this coupling significantly simplifies the calculations of primordial abundances and 
observable signals at colliders and underground detectors, making the resulting predic- 
tions more certain. The primordial abundance is governed mainly by the annihilation 
of S-particles via s-channel Higgs exchange. If the mass of the Higgs is not too close 
to 2m s , the observed abundance of dark matter is achieved by a most natural choice 
for the coupling, A ~ 0(0.1 — 1). 

These large values of A lead to nuclear elastic cross sections which are ~ 10~ 43 
cm 2 in size. This is slightly below the limits of sensitivity of the DAMA and CDMS 
experiments, but is likely to be detectable at future experimental facilities. Thus, the 
direct detection of S particles is not yet ruled out, and is easily feasible within the near 
future! Our minimal model is unable to reconcile the DAMA and CDMS results, and 
so predicts that one or the other is incorrect. 

Cross sections this large also lead to a significant flux of high-energy neutrinos 
generated by S annihilation at the centers of the earth and the sun, at levels potentially 
accessible to large neutrino detectors. We give here only an estimate of the expected 
flux, and we believe the encouraging results motivate further work to simulate the 
expected neutrino and muon spectra and intensities. It would also be intriguing to 
repeat the numerical analysis of ref. JTO , updated with the Higgs masses favored by 
modern collider results. 

Large values for A may also lead to the significant missing-energy signals at colliders, 
corresponding to the production of a pair of S'-particles. This signal is unlikely to be 
seen at hadronic machines, where the background events (two jets plus missing energy) 
will have much larger cross-sections. Lepton colliders, such as LEP or NLC are needed. 
The possible existence of S-particles with a mass smaller than half of the Higgs boson 
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mass poses a significant threat to Higgs searches at the Tevatron and LHC. Indeed, 
in this case the invisible decay of the Higgs boson can be many orders of magnitude 
larger than the search modes, making the Higgs effectively invisible. For intermediate 
Higgs masses, up to 140 GeV and higher, the pair production of W bosons introduces a 
nonnegligible visible width and so saves the Higgs searches. It is important that in the 
dangerous domain of parameter space, 2m s < < 140 GeV, the size of the elastic 
cross-section, a e \, is limited from below to be larger than 10 -44 cm 2 . If this were indeed 
the case, the next generation of the underground detectors would likely discover the 
recoil signal. 

Because the self-coupling, A s , is unconstrained by the abundance condition, one 
might hope to use models of this sort to produce a strongly-interacting dark matter 



candidate, a possibility recently advocated in ref. [[12], |I4"|| . Unfortunately, this proposal 
would also require a rather small mass for the S-scalars (1 GeV or less), and since A 
cannot be small for such masses, one is led to a fine-tuned (ppm) cancellation between 
\Vg W and the "bare" mass m 2 ,. Worse, having m s be 500 MeV or smaller pushes A to 
the very limits of our perturbative analysis. 

It remains an interesting challenge to construct a non-minimal model of strongly- 
interacting dark matter along the present lines. One possibility is to suppress A to 
evade the fine-tuning problem, and set the freeze-out abundance by a set of non- 
renormalizable operators of dimension 6 and higher. Another possibility for having 
light and strongly interacting scalar dark matter requires A to be very small, scalar 
particles to be out of equilibrium at all temperatures, and their abundance to be fine- 
tuned (for example, using inflationary scenarios). The suppression of the coupling to 
the Higgs, strong self-interaction and a low mass for 5* could be simultaneously achieved 



in models where 5* is composite plf , in which case there is unlikely to be significant 
impact on Higgs physics. 

Our pursual here of a more model-independent approach to the physics of dark 
matter particles is ultimately stimulated by the maturity of current dark-matter de- 
tection efforts. As the data erodes the large parameter space of the popular neutralino 
models, or if a dark matter candidate is eventually found, a reliable interpretation of 
the experiments can only be found by comparing with a wide variety of explanatory 
models. Our proposal in this paper marks the absolute minimum of complexity re- 
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quired by such a model, and from this minimality springs its predictiveness. It is also 
generic in the sense that it is the most general renormalizable model consistent with 
the assumed particle content and symmetries. As such, it captures the implications of 
any models whose implications for dark matter arise from a low energy spectrum which 
contains only these particles and symmetries. We believe that further model-building 
efforts are in order as experimental developments progress. 



Acknowledgements 

We thank G. Couture, T. Falk, L. Kofman, G.Mahlon, K. Olive, S. Rudaz and 
M. Voloshin for usefull discussions and comments. M.P. wishes to thank the Physics 
Departments of McGill University and the University of Quebec at Montreal, where 
a significant portion of this work was carried out. M.P. and T. t V.'s research was 
funded by the US Department of Energy under grant number DE-FG-02-94-ER-40823, 
and C.B. acknowledges the support of NSERC (Canada), FCAR (Quebec) and the 
Ambrose Monell Foundation. 



23 



References 



[1] M. S. Turner, Phys. Rept. 333-334 (2000) 619. 

[2] See e.g. G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 (1996) 
195, L. Bergstrom, Rept. Prog. Phys. 63 (2000) 793. 

[3] H. Pagels and J. R. Primack, Phys. Rev. Lett. 48 (1982) 223; 

J. Ellis, J.S. Hagelin, D.V. Nanopoulos, K.A. Olive and M. Srednicki, Nucl. Phys. 
B238 (1984) 453. 

For a recent study, see: J. Ellis, T. Falk, G. Ganis and K. A. Olive, Phys. Rev. 
D62 (2000) 075010. 

[4] D. Caldwell, Nucl. Phys. B70 (Proc. Suppl.) (1999) 43. 

[5] R. Bernabei et al, preprint ROM2F/2000/01 and INFN/AE-00/01; 
R. Bernabei et al, Phys. Lett. B450 (1999) 448. 

[6] R. Abusaidi et al, Phys. Rev. Lett. 84 (2000) 5699. 

[7] M. Pospelov and T. ter Veldhuis, Phys. Lett. B480 (2000) 181. 

[8] M. Veltman and F. Yndurain, Nucl. Phys. B325 (1989) I. 

[9] V. Silveira and A. Zee, Phys. Lett. B161 (1985) 136. 

[10] J. McDonald, Phys. Rev. D50 (1994) 3637. 

[11] R. Flores and J. M. Primack, Ap. J. 427 (1994) LI; 
B. Moore, Nature 370 (1994) 620; 

B. Moore et a/., M.N.R.A.S. 310 (1999) 1147, Ap. J. 524 (1999) L19; 
A. A. Klypin et a/., Ap. J. 522 (1999) 82, Ap. J. 554 (2001) 903. 

[12] N. Spergel and P.J. Steinhardt, Phys. Rev. Lett. 84 (2000) 3760. 

[13] N. Yoshida et a/., Ap. J. 544 (2000) L87. 

[14] M.C. Bento et.al, Phys. Rev. D62 (2000) 041302. 



24 



[15] R. S. Chivukula and M. Golden, Phys. Lett. B267, (1991) 233; 

R. S. Chivukula, M. Golden and M. V. Ramana, Phys. Lett. B293, (1992) 400; 

T. Binoth and J.J. van der Bij, [hep-ph/9511467| ; 

T. Binoth and J. J. van der Bij, Z. Phys. C75 (1997) 17 

R. Akhoury, J. J. van der Bij and H. Wang, |hep-ph/0010187 . 

[16] S. Coleman, Nucl. Phys. B262, (1985) 263. 

For a recent discussion, see: D. A. Demir, Phys. Lett. B450, 215 (1999). 



[17 
[18 
[19 
[20 
[21 
[22 
[23 

[24; 

[25 
[26 
[27 
[28 
[29 
[30 
[31 
[32 



L. Kofman, A. Linde and A. A. Starobinsky, Phys. Rev. D56 (97) 3258. 
E.W. Kolb and M.S. Turner, The Early Universe, Addison- Wesley, 1990. 
K. Griest and D. Seckel, Phys. Rev. D43 (1991) 3191. 
M. Voloshin, Sov. J. Nucl. Phys. 44 (1986) 478. 
A. Faraggi and M. Pospelov, ^iep-ph/0008223 . 



R. Bernabei et al, Phys. Lett. B389 (1996) 757. 

PF. Smith et al, Phys. Lett. B376 (1996) 299. 

M. Beck et al, Phys. Lett. B336 (1994) 141. 

A. Morales et al, Phys. Lett. B489 (2000) 268. 

M. Mori et al, Phys. Rev. D48 (1993) 5505. 

M. Ambrosio et al, Phys. Rev. D60 (1999) 082002. 

M.M. Boliev, et al, Nucl. Phys. B48 (Proc. Suppl.) (1996) 83. 

E. Andres et al, Nucl. Phys. Proc. Suppl. 91 (2000) 4232. 

A. Okada, Nucl. Phys. Proc. Suppl. 91 (2000) 423. 

M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Phys. Lett. B78 (1978) 264. 
A Zhitnitsky, Phys. Rev. D55 (1997) 3006. 



25 



[33] A. E. Faraggi, K. A. Olive and M. Pospelov, Astropart. Phys. 13 (2000) 31. 

[34] J. Silk, K.A. Olive, and M. Srednicki, Phys. Rev. Lett. 55 (1985) 257; 
M. Srednicki, K.A. Olive, and J. Silk, Nucl. Phys. B279 (1987) 804. 

[35] S. Ritz and D. Seckel, Nucl. Phys. B304 (1988) 877. 

[36] J. Edsjo, Nucl. Phys. Proc. Suppl. 43 (1995) 265. 

[37] T.K. Gaisser and T. Stanev, Phys. Rev. D31 (1985) 2770. 

[38] M. Honda and M. Mori, Prog. Theor. Phys. 78 (1987) 963. 

[39] C. Burgess, J. Matias and M. Pospelov, |hep-ph/9912459| . 

[40] H. Haber, D. Dicus, M. Drees, and X. Tata, Phys. Rev. D37 (1987) 1367; 
K. Griest and H. E. Haber, Phys. Rev. D37 (1988) 719; 
A. Djouadi, J. Kalinowski, and P.M. Zerwas, Z. Phys. C57 (1993) 569; 
J. L. Lopez, D. V. Nanopoulos, H. Pois, X. Wang, and A. Zichichi, Phys. Rev. 
D48 (1993) 4062. 

[41] J. F. Gunion, Phys. Rev. Lett. 72 (1994) 199. 

[42] S. P. Martin and J. D. Wells, Phys. Rev. D60 (1999) 035006. 

[43] O. J. Eboli and D. Zeppenfeld, Phys. Lett. B495 (2000) 147. 

[44] P. I. Kemenes, talk given at the XXXth International Conference on High Energy 
Physics, Osaka 2000; Searches for Higgs Bosons: Preliminary Combined Results 
using LEP Data Collected at Energies up to 202 GeV, the ALEPH, DELPHI, 
L3 and OPAL Collaborations (P. Bock et al), preprint CERN-EP-2000-055, Apr 
2000; 

This limit has recently changed appreciably, since the bound from previous running 
gave a limit of 94.5 GeV (ALEPH). See: 

ALEPH Collaboration (R. Barate et al.), Phys. Lett. B465 (1999) 349; 
DELPHI Collaboration (P. Abreu et al.), Phys. Lett. B459 (1999) 367; 
L3 Collaboration (M. Acciarri et al.), Phys. Lett. B485 (2000) 85; 
OPAL Collaboration (G. Abbiendi et al.), Eur. Phys. J. C12 (2000) 567. 



26 



